Phase Diagram of the Two-Channel Kondo Lattice 
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The phase diagram of the two-channel Kondo lattice model is examined with a Quantum Monte 
Carlo simulation in the limit of infinite dimensions. Commensurate (and incommensurate) antiferro- 
magnetic and superconducting states are found. The antiferromagnetic transition is very weak and 
continuous; whereas the superconducting transition is discontinuous to an odd-frequency channel- 
singlet and spin-singlet pairing state. 
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Introduction. A number of Heavy Fermion materials 
display highly unusual superconductivity (HFSCs) It 
seems unlikely that conventional superconductivity co- 
exists with the strong local Coulomb correlations neces- 
sary to enhance the electronic mass. Indeed, the spe- 
cific heat jump at the transition scales with the normal 
state specific heat making it clear pairing is between the 
heavy electrons (with order 100-1000-fold mass enhance- 
ment). The strong coulomb correlations present no prob- 
lem for unconventional superconducting order parame- 
ters with either spatial or temporal nodes (so called 
"odd-frequency pairing) [gpQ . Such an interpretation is 
supported by two sets of data: (i) power laws observed 
in physical properties below the superconducting tran- 
sition which contrast with the activated behavior 
of the conventional (nodeless) s-wave order of, e.g., alu- 
minum; (ii) the complex superconducting phase diagrams 
of UPt 3 and Ui^ThrBeis (and possibly UBei 3 itself, for 
which the penetration depth displays evidence of a sec- 
ondary transition) 0J^Q ■ The superconductivity can co- 
exist with antiferromagnetism, which is usually commen- 
surate, as in UPt 3 , URu 2 Si 2 , UPd 2 Al 3 , and UNi 2 Al 3 g), 
or compete with it as in CeCu 2 Si 2 |9j . Finally, at least 
in UBei 3 |l(J and CeCu 2 Si 2 P], the superconductivity 
arises in a normal state which is clearly not described as 
a Fermi liquid. In each of these materials above T c the 
linear specific heat rises with decreasing temperature, the 
resistivity is approximately linear in T, and the residual 
resistivity at T c is high (typically 80-100 /i— Q,-cm in the 
best samples of UBei 3 ). 

In this paper, we provide the first calculations of the 
phase diagram for the two-channel Kondo lattice in in- 
finite spatial dimensions. We find second order antifer- 
romagnetic and first order odd frequency superconduct- 
ing phase transitions. Coexistence of commensurate (and 
incommensurate) antiferromagnetism with superconduc- 
tivity is in principle possible. We present several pos- 
sible routes to account for the multiple superconduct- 
ing phases observed in real materials. Taken together 
with the non-Fcrmi liquid (NFL) paramagnetic phase dis- 
cussed in previous publications and earlier sugges- 
tions that the two-channel lattice may describe UBei 3 



and other heavy fermion systems 0, our work estab- 
lishes the two-channel Kondo lattice model as possessing 
the key ingredients needed to explain heavy fermion su- 
perconductivity. 

Motivation A possible model for some heavy fermion 
compounds is the two-channel Kondo lattice which con- 
sists of two identical species of non-interacting elec- 
trons antiferromagnetically coupled to an adjacent lat- 
tice of spin 1 /2 Kondo moments ^ . This model displays 
Non-Fermi liquid behavior because of the overcompen- 
sation of the Kondo spins by the conduction electrons, 
first pointed out in the single site two-channel model 
by Nozieres and Blandin JlT| where extended screening 
clouds form which retain a spin 1/2 character. As a re- 
sult, the paramagnetic state is degenerate and the exci- 
tation spectrum non-Fermi liquid like. This overscreen- 
ing results in an interchannel pairing mechanism Q in 
the two-channel impurity problem which tends to favor 
spin-singlet channel-singlet odd frequency superconduct- 
ing fluctuations For the lattice, we find that this 
overscreening generates novel antiferromagnetic superex- 
change between the Kondo spins which, taken together 
with the RKKY exchange, favors antiferromagnetism in 
the lattice model close to half filling of the bands. The na- 
ture of the resulting superconducting and magnetic tran- 
sitions, as well as the competition between them, is the 
subject of this manuscript. . 

Model The Hamiltonian for the two-channel Kondo 
lattice is 



H 



t* 



2Vd 



+ h.c. 



<ij >,a,er 



(1) 



where c\ a (ci ya ^) creates (destroys) an electron on site 
i in channel a = 1, 2 of spin a, Si is the Kondo spin on 
site i, and Si^ a are the conduction electron spin operators 
for site i and channel a. The sites i form an infinite- 
dimensional hypercubic lattice. Hopping is limited to 
nearest neighbors with hopping integral t = t* /2\fd\ the 
scaled hopping integral t* determines the energy unit and 
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is set equal to one (t* = 1). Thus, on each site the 
Kondo spin mediates an interaction between the two dif- 
ferent channels. This problem is non-trivial, and for the 
region of interest in which J > and T <C J, t* it is de- 
scribable only with non-perturbative approaches. Clearly 
some simplifying method which allows for a solution of 
the lattice problem in a non-trivial limit is necessary. 

Formalism and Simulation. Such a method was pro- 
posed by Metzner and Vollhardt |r| who observed that 
the renormalizations due to local two-particle interac- 
tions become purely local as the coordination number 
of the lattice increases. A consequence is that the so- 
lution of most standard lattice models may be mapped 
onto the solution of a local correlated system coupled to 
an effective bath that is self-consistently determined. We 
refer the reader to recent reviews for further details and 
references on the method [jl4| . 

In order to solve the remaining impurity problem, we 
use the Kondo impurity algorithm of Fye and Hirsch 
fl5|| , modified to simulate the two-channel problem [fl6|| . 
We simulated the model for a variety of fillings and 
values of the magnetic exchange J (0 < N < 1 for 
J = 0.75,0.625,0.5,0.4). Error bars on measured quan- 
tities are less than 6% for the results presented here. A 
sign problem was also encountered in the QMC process. 
It becomes severe for large J and also for low fillings since 
lower temperatures are required to access the physically 
interesting regime. This precludes a systematic study of 
the model with J > 0.8 and filling N < 0.5. 

Both one and two-particle properties arc natural prod- 
ucts of the QMC simulation. To calculate the spin sus- 
ceptibility, it was necessary to measure the local three- 
by-three matrix of the susceptibility including both the 
Kondo spin and conduction band spin fluctuations. This 
may then be inverted to calculate the associated irre- 
ducible vertex function and the corresponding lattice sus- 
ceptibility in the usual way 

The situation for the superconductivity is a bit more 
complicated. Since there are no charge fluctuations in- 
volving the Kondo spin degree of freedom, only the two 
conduction channels contribute to the pair-field suscep- 
tibility. We can then look for pairing instabilities in sin- 
glet and triplet channels for both spin and channel spin. 
Motivated by the result for the impurity model [^|^] we 
have restricted our attention to interchannel and inter- 
spin particle-particle diagrams as shown in Fig. 1. It is 
possible to make two independent combinations of these 
diagrams x±(ia; n ,£w m ) = Xn(^Ww m ) ±xi2(i^ n ,iw m ), 
from which we construct a quartet of spin and channel, 
singlet and triplet pair-field susceptibilities given by 

PssCs = T /_ (iu n )x- (iw n , iu m )f- (iu> m ) (2) 
nm 

Pstct = f+(ivn)x+(ivn,iv m )f+(iu m ) (3) 

nm 



Pstcs = T^x-fan, ium) (4) 

n m 

Pssct =T^2x+( iuj n,iu m ) ■ (5) 

nm 

Here f±(iui n ) are °dd functions of Matsubara frequency 
used to project out the odd-frequency pairing, and, for 
example, PssCs{T) is the spin-singlet channel-singlet 
pair-field susceptibility. 

To determine the form of f±(iu) n ), we employ the pair- 
ing matrix formalism of Owen and Scalapino Jl8| . Here 
we represent each of x± m a two-particle Dyson equa- 
tion and extract the irreducible vertex functions T± . The 
pairing matrices are then 

M ± (iu n ,iuj m ) = ^Jx± ( i V n ) r± (i LU n , «w m ) y/ x± {i^m ) ■ 

(6) 

where x±(iw n ) are the particle-particle diagrams in Fig. 1 
without vertex corrections. f±(iu> n ) is the eigenvector 
corresponding to the dominant eigenvalue of M± (that 
with the largest absolute value), 

Results. In this model antiferromagnetism is driven 
by both RKKY interactions and a novel type of superex- 
change. The latter arises from hopping between adja- 
cent spin 1/2 screening clouds, whose overall spin is de- 
termined by the conduction electrons; the Pauli princi- 
ple forbids hopping unless neighboring spins in the same 
channel are antiparallel. As a result, for large J, the 
superexchange goes as ~ (t*) 2 /J. For conduction band 
fillings close to N — 1 both the RKKY (evaluated at 
nearest neighbor sites) and the superexchange favor anti- 
ferromagnetism (the RKKY exchange remains antiferro- 
magnetic until N < 0.5), and an antiferromagnetic tran- 
sition results, as shown in Fig. 2. Due to screening of 
the local moments by the conduction spin, the transi- 
tion is very weak, as measured by the full susceptibil- 
ity. Specifically, xaf is not significantly enhanced over 
the bulk susceptibility xf until T > T/v- However, the 
f-electron contribution to the susceptibility shows a pro- 
tracted scaling region. Here screening effects non-linear 
feedback that reduces the susceptibility exponent 7 from 
the mean field value 7 = 1. 7 increases with doping 
(N < 1), and the transition becomes incommensurate as 

To explore superconductivity, it is necessary to find the 
frequency form factors f± , which, as discussed above, are 
the eigenvectors corresponding to the dominant eigenval- 
ues of the pairing matrices M±. As the temperature is 
lowered, the dominant eigenvalue first becomes large (di- 
vergent) and negative, and then abruptly switches to a 
large and positive value at the transition. This happens 
first in M_, and the corresponding eigenvector of M_ 
is plotted in the inset to Fig. 3. It can be fit quite ac- 
curately to the form T /2uj n as shown by the solid line, 
which corresponds to the form factor of Rcf. H . Thus, we 
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use f-(iuj n ) — T/2u n to project out the odd-frequency 
pair-field susceptibilities shown in Figs. 3 and 4. (Other 
form factors f(iuj n ) = tanh(To/a; n ) and sign(oj n ) were 
tried and produce qualitatively similar results.) 

The first transition is found in the spin-singlet channel- 
singlet pairing combination, as shown in Fig. 3. (Note 
that this pair state is of necessity even in parity, so that 
the odd-frequency condition is required to satisfy the 
Pauli principle, in contrast with the spin-parity selec- 
tion rules in the single channel case |j|,f|].) To interpret 
this result, remember that the inverse pair field suscepti- 
bility is proportional to the curvature of the free energy 
f(AssCs) as a function of the pairing order parameter 
1/PssCs oc d 2 f(A SsCs )/dA 2 SsCs . Thus, if P SsC s < a 
thermodynamic instability of the system is present. The 
associated transition cannot be continuous, since this re- 
quires the free energy in order parameter space to be- 
come fiat (ie. PssCs diverges) so that the order parame- 
ter may change continuously. Thus, we identify this as a 
discontinuous transition. Furthermore, if PssCs{T*) = 0, 
then T* is a lower bound to the transition temperature 
since when Ps s Cs{T — T* = Q~) = 0~ the curvature of 
f(AssCs) is divergent and negative. I.e. the free en- 
ergy displays a downward cusp which would compel the 
order parameter and the free energy to change discon- 
tinuously at T*, which involves an infinite energy at the 
transition. Hence, the actual transition occurs at a tem- 
perature T c > T*. 

Several remarks are in order about this unusual super- 
conductivity: (a) Fig. 3 suggests that the transition is 
degenerate at the zone center and the zone corner, and in 
fact PssCs vanishes simultaneously over the whole zone, 
and hence this is a locally driven transition! This degen- 
eracy may be lifted when the order parameter becomes 
finite, otherwise, it will be lifted in finite dimensionality 
as the superfluid stiffness T s (analogous to spin stiffness 
in a magnet) must vanish in infinite d. For finite d with 
T s > 0(< 0), a local free energy minimum will be found 
at q = Q(q ^ 0). (b) The vanishing of the local pair sus- 
ceptibility here at T* contrasts with the impurity model 
in which it diverges logarithmically as T — > [j|rc|. (c) 
The excitation spectrum of such a transition may well be 
highly exotic 0). 

As shown in Fig. 4, the ground state of the system 
may be superconducting or magnetic which may coex- 
ist or compete. Our work gives only estimates of the 
transition temperatures; detailed exploration within the 
ordered phases will answer this question definitively. In 
general superconductivity will occur first (at the highest 
transition temperature) away from half filling and anti- 
ferromagnetism will occur first near half filling. However, 
for values of J > 0.75, we find that superconductivity 
occurs first even at half filling (the minus sign problem 
precluded a systematic study for these values of J so they 
are not presented in Fig. 4). We also find that, T* cx To, 
with the constant of proportionality sa 0.5 weakly depen- 



dent on J; whereas T/v appears to depend upon both J 
and To. 

Speculation and Interpretation We offer two compati- 
ble interpretations of the result T* w 0.5Xb- First, consis- 
tent with the fact that the transition is locally driven, the 
transition temperature is proportional to the Kondo scale 
To, the temperature characteristic of the local dynamics. 
Second, T rs 0.5To is also the temperature at which the 
slope of the real part of the self energy becomes one , 
so that the quasiparticle renormalization factor diverges. 
This result suggests that when the system cannot form a 
Fermi liquid due to a large residual scattering rate (which 
also yields a large residual entropy) , it forms a supercon- 
ducting state to quench the entropy. The antiferromag- 
netic transition temperature, on the other hand, depends 
upon both J (through the intersite exchange) and upon 
To (through moment screening and superexchange) . As 
noted above, we cannot determine whether these states 
coexist without a detailed exploration within the ordered 
phases. At half filling, where the antiferromagnetism pro- 
duces an insulating phase, clearly the superconductivity 
will be suppressed for J < 0.75. 

Our results offer a number of routes to be explored for 
explaining the complex superconducting phase diagrams 
of UPt3 and Ui-^Th^Beia: 1) Competition between 
phases with a multi-point irreducible star. As men- 
tioned above, for T s < in finite d, q ^ pairs are fa- 
vored. In this simple case of a bipartite lattice, the zone 
corner has a one-point irreducible star. For lattices with 
frustration, such as hexagonal UPt3 and face-centered 
cubic UBei3, it is possible to produce multi-point irre- 
ducible stars for staggered order parameters which can 
then have multiple phases |L9]]. 2) Competition Be- 
tween Different q values. As noted above, T*, the 
lower bound for the first order superconducting transi- 
tion temperature is independent of q in our calculations. 
Thus multiple phases may thus correspond to supercon- 
ducting transitions with different q values. 3)Possible 
Instability of Spin Triplet/Channel Triplet Pair- 
ing. At yet lower temperatures than those identified in 
Fig. 4, we observe a sign change in the pair field suscep- 
tibility associated with spin triplet-channel triplet odd- 
frequency pairing. Hence it is possible that the compe- 
tition between this triplet-triplet and the singlet-singlet 
odd frequency pairing may explain the complex phase 
diagrams. 

In closing, we note that this superconducting transi- 
tion can agree with experiment only if it is weakly first 
order; this is plausible given the rapid change in free en- 
ergy curvature at the order parameter origin. Detailed 
investigations in the ordered phase will resolve this issue 
and whether any of the above scenarios can describe the 
complex phase diagrams of UPt3 and Ui-xThaBeia. 
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FIG. 1. Particle-particle interchannel opposite spin dia- 
grams which contribute to the pair-field susceptibility. Here 
1 and 2 label the channel and j and J. the spin. 
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FIG. 2. Full antiferromagnetic xaf and ferromagnetic 
Xf susceptibilities of the two-channel Kondo lattice with 
J = 0.75 and N = 1.0. xaf shows very little enhancement 
until T ~ Tjv. As shown in the inset, the Kondo spin con- 
tribution to the susceptibility Xaf(T) displays a protracted 
scaling region. The solid line is a fit of 1/Xaf( t ) t0 the form 
a(T — Tjv) 7 . The large difference between xaf and Xaf(T) as 
well as the reduced exponent 7 < 1 indicate the importance 
of screening of the local spin in this transition. 
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FIG. 3. Odd-frequency (spin-singlet, channel-singlet) 
pair-field susceptibility at both the zone center and corner. 
At low temperatures T = T* < T c PssCs{T) becomes nega- 
tive, indicating a thermodynamic instability to the formation 
of a pairing state. The lack of a divergence of PssCs(T) rules 
out a continuous transition, hence the system must undergo 
a discontinuous transition. In the inset, the dominant eigen- 
vector of the pairing matrix M_ is plotted versus Matsubara 
frequency. It can be fit to the form T/2u> n (solid line), which 
corresponds to the frequency variation of the gap function 
found in Ref.4 



FIG. 4. Phase diagrams of the two-channel Kondo lattice 
for various values of J. The solid lines are fits to the data. 
For each value of J shown the antiferromagnetic transition 
occurs first (at a higher temperature) at half filling N = 1 
(albeit very weakly when J = 0.75). Tm falls quickly with 
doping N < 1, with the transition becoming incommensurate 
near Tjv — ► 0. Away from half filling, superconductivity oc- 
curs first above a temperature T* > T c is roughly equal to 
CT , (C = 0.43, 0.51, and 0.58, for J = 0.5,0.625,0.75 respec- 
tively) consistent with the local nature of the transition. T* 
should be regarded as a lower bound on the superconducting 
transition temperature T c (see text). 
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